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Single fluorescent molecules in a flowing 
sample stream (14) are distinguished and identified 
using only a single laser excitation wavelength 
(38). A sample stream (14) is formed containing 
a dilute mixture of single molecule fluorophores 
(12), wherein each one of the fluorophores is 
serially ordered in the sample stream (14). The 
sample stream (14) is illuminated with a single 
excitation wavelength laser (38) effective to excite 
each fluorophore one at a time. Fluorescence 
emission photons from each said fluorophore are 
detected (48). A burst size is determined for 
each fluorophore to identify each fluorophore. A 
pulsed laser may be used (38), where burst size 
and an intra-burst fluorescence decay rate for each 
fluorophore are determined simultaneously from 
the detected fluorescence emission photons. The 
burst size and the decay rate are correlated to 
identify each fluorophore. 
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SINGLE MOLECULE IDENTIFICATION USING 
SELECTED FLUORESCENCE CHARACTERISTICS 



RELATED CASES 

This application claims the benefit of the priority date of provisional 
application S.N. 60/065,365, filed November 12, 1997. 

BACKGROUND OF THE INVENTION 

This invention relates to flow cytometry and, more particularly, to 
applications of flow cytometry using single molecule identification. This 
invention was made with government support under Contract No. W-7405- 
ENG-36 awarded by the U.S. Department of Energy. The government has 
certain rights in the invention. 

Light-induced fluorescence detection of single molecules in liquid 
solution was first accomplished several years ago, and applications for single 
molecule detection (SMD) in the analytical, environmental, and biomedical 
sciences are beginning to emerge. Overviews of single molecule detection in 
solution are presented in, e.g., R.A. Keller et al., 5Q Appl. Spectrosc, pp. A12- 
A32 (1996), and P.M. Goodwin et al., 2& Accounts Chem. Res., pp. 607-613 
(1996). In general, an operable system will include a dilute stream of 
separated individual molecules with known fluorescent characteristics that are 
excited one at a time by a light source, where the resulting emitted photons 
are detected. One approach to single molecule detection used herein is 
based on flow cytometry, where the analyte solution is delivered into a rapidly 
flowing sheath fluid and hydrodynamically focused into a narrow sample 
stream. Yet another approach is based on analyte movement through a 
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capillary, such as described in U.S. Patent 5,209,834, issued May 11, 1993, to 
Shera, and incorporated herein by reference. 

In flow cytometry, a sample stream passes through the center of a 
probe volume, defined by the diameter of the focused excitation laser beam 
5 and a spatial filter placed in the image plane of a light collecting objective. 

Single fluorescent molecules are detected by the bursts of photons that they 
emit as they flow through the detection volume one-at-a-time. Hydrodynamic 
focusing of the sample stream by the sheath fluid ensures that the entire 
sample stream passes through the center of the excitation laser, so that single 
10 molecules delivered into the flow cell are detected with an efficiency 

exceeding 90%. See, e.g., P.M. Goodwin et al., "Progress toward DNA 
sequencing at the single molecule level," 41 Exp. Tech. Phys., pp. 279-294 
(1995), incorporated herein by reference. 

Some of the applications under development for efficient single 
15 molecule detection in flow include DNA fragment sizing, DNA sequencing, 
counting and sorting of single molecules, and detection of probe-target 
binding. A number of applications for SMD in solution require one to 
distinguish between different fluorophores present in a mixture. For example, 
in one approach to DNA sequencing (U.S. Patent 4,962,037, issued October 
20 9, 1990, and incorporated herein by reference), each base is labeled with a 
different fluorescent probe, and a rapid, efficient method is needed to identify 
these fluorophores. In yet another approach, only two fluorescent probes are 
required (U.S. Patent 5,405,747, issued April 11. 1995, and incorporated 
herein by reference) to reduce the number of distinguishing characteristics 
25 that are required to be identified. 

Several techniques have been developed to distinguish between 
different single molecules in solution. One technique employs two or more 
detection channels to identify single molecules in a multi-component mixture, 
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based upon differences in excitation and emission wavelengths of the 
fluorophore labels. See, e.g., Soper et al. f "Detection and identification of 
single molecules in solution;* 2 J. Opt Soc. Am. B, pp. 1761-1769 (October 
1992), incorporated herein by reference. Typically, the detection volume is 

5 probed by collinear laser beams at different wavelengths, corresponding to 
the excitation maxima of each of the fluorophores. A beam splitter directs the 
fluorescence from each fluorophore to a separate filter/detector channel. This 
method requires that the emission bands of the fluorophores be sufficiently 
separated to minimize cross-talk between the respective channels. 

10 Furthermore, separating the fluorescence signal into distinct detection 

channels increases the complexity of the instrumentation and can reduce the 
overall detection efficiency. 

An alternative approach that requires only a single detection channel 
and is applicable to molecules with similar spectroscopic properties exploits 

15 the difference in fluorescence lifetimes of the fluorophores. Single Rhodamine 
6G (R6G) molecules in flow have been distinguished from Rhodamine B (RB) 
molecules and from tetramethyl rhodamine isothiocyanate (TRITC) molecules 
by applying a maximum likelihood function to time-correlated single-photon 
counting (TCSPC) fluorescence decay measurements of individual bursts in a 

20 mixed sample (Zander et al.,63 Appl. Phys. B, pp. 517-523 (1996); J. 
Enderlein et al. v 2ZQ Chem. Phys. Lett,, pp. 464^70 (1997)). Also, a 
rhodamine derivative JA169 has been distinguished from a carbocyanine dye 
Cy5 using this technique (M. Sauer et at., 65 Appl. Phvs. B. No. 3 . pp. 427- 
431 (August 1997). 

25 Objects, advantages and novel features of the invention will be set forth 

in part in the description which follows, and in part will become apparent to 
those skilled in the art upon examination of the following or may be learned by 
practice of the invention. The objects and advantages of the invention may be 



SUBSTITUTE SHEET (RULE 26) 



WO 99/24797 PCT/US98/241 1 8 

4 

realized and attained by means of the instrumentalities and combinations 
particularly pointed out in the appended claims. 

SUMMARY OF THE INVENTION 
5 To achieve the foregoing and other objects, and in accordance with the 

purposes of the present invention, as embodied and broadly described herein, 
this invention comprises a method for identifying single fluorescent molecules 
in a flowing sample stream. A sample stream is formed containing a dilute 
mixture of single molecule fluorophores, wherein each one of the fluorophores 
10 is serially ordered in the sample stream. The sample stream is illuminated 
with a single excitation wavelength laser beam effective to excite each 
fluorophore. Emitted photons from each said fluorophore are detected. In 
one embodiment, a burst size is determined for each fluorophore to identify 
each fluorophore. In another embodiment, the laser beam is formed of pulses 
15 and a burst size and an intra-burst fluorescence decay rate (/) for each 

fluorophore are determined from fluorescence emission photons. The burst 
size and the decay rate are correlated to identify each fluorophore. 

BRIEF DESCRIPTION OF THE DRAWINGS 
20 The accompanying drawings, which are incorporated in and form a part 

of the specification, illustrate the embodiments of the present invention and, 
together with the description, serve to explain the principles of the invention. 
In the drawings: 

FIGURE 1 schematically depicts apparatus for identifying single 
25 fluorescent molecules. 

FIGURES 2A-2D graphically depict single molecule fluorescence 
photon burst data from a photolysed ultrapure water blank, a dilute sample 
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stream of TRITC t a dilute sample stream of R6G, and a dilute sample stream 
with a mixture of TRITC and R6G, respectively. 

FIGURES 3A and 3B graphically depict single molecule fluorescence 
burst duration distributions measured from dilute sample streams of TRITC 
and R6G, respectively. 

FIGURES 4A-4F graphically depict fluorescence photon burst size 
distributions (BSDs) compiled from data collected from a dilute sample stream 
of R6G, of TRITC and a mixture of R6G and TRITC (FIGURES 4A, 4C, 4E) 
and corresponding BSDs that are time filtered (FIGURES 4B, 4D, 4F). 

FIGURES 5A and 5B are correlated single molecule burst size and 
intra-burst fluorescence decay rate (/") measurements for a sample stream 
containing approximately equal amounts of TRITC and R6G and a sample 
stream containing mostly R6G. 

DETAILED DESCRIPTION OF THE INVENTION 
In accordance with the present invention, single molecule identification 
(SMI) is based on differences in detected fluorescence burst intensities from 
fluorescent molecules. These intensity differences are a consequence of the 
different photophysical properties (absorption cross-sections, fluorescence 
quantum yields, and fluorescence emission spectra) of different fluorescent 
molecules. Unlike discrimination methods based on separate excitation 
wavelengths and emission detection channels or by fluorescence decay rates, 
burst intensity discrimination requires that the molecules experience a similar 
excitation irradiance during their transit across an excitation laser beam. The 
excitation laser beam may be either a pulsed laser beam or a continuous- 
wave laser beam. To accomplish this, the molecules are constrained to a 
sample stream diameter within the focused excitation laser. As with 
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fluorescence decay rate measurements, this technique can be applied to 
molecules with similar spectroscopic properties and requires only a single 
excitation wavelength and fluorescence emission detection channel. In 
addition, a multiplex technique, wherein individual fluorescent molecules are 
distinguished by simultaneous measurement of fluorescence burst intensity 
and intra-burst fluorescence decay rate, provides a relatively high probability 
of individual fluorescent molecule identification. When intra-burst 
fluorescence decay rate is determined, the excitation laser beam must be a 
pulsed laser beam. 

EXPERIMENTAL SECTION 
SMD in a flowing sample stream has been detailed in several 
publications, as noted in the Background of the Invention. The experimental 
techniques and results presented herein are also reported in Van Orden et ah, 
"Single-Molecule Identification in Flowing Sample Streams by Fluorescence 
Burst Size and Intraburst Fluorescence Decay Rate," J£LAnai Chem. No. 7, 
pp. 1444-1451 (April 1, 1998), incorporated herein by reference. An operable 
system requires only a sample stream having an ordered flow of molecules 
that pass through the center of an excitation and detection volume. The 
experimental apparatus described herein provides increased efficiency for the 
detection of single molecules, and is a preferred apparatus, but the present 
invention is not limited to flow cytometry. Other devices are available that 
provided an ordered stream of molecules, e.g., capillary flow and 
microchannels, are known to persons of ordinary skill in the detection of single 
molecules and are within the scope of the present invention. 

To obtain the results reported below, apparatus similar to that 
schematically shown in Figure 1 was used. The parameters specified below 
are exemplary only and not intended to limit the scope of the present 
invention. Solutions 12 of, e.g., R6G and TRITC fluorophore molecules 



SUBSTITUTE SHEET (RULE 26) 



WO 99/24797 



PCIYUS98/24118 



7 

dissolved in ultrapure water, were delivered 14 into a sheath fluid IS, flowing 
through a square bore flow cell 18 (250 x 250 pm 2 ), from fused silica capillary 
22 (10-pm i.d., 140-pm o.d.). The outlet tip 24 of the capillary was etched to a 
narrow taper using hydrofluoric acid to facilitate smooth laminar flow of the 
sheath fluid around the capillary. The inlet end 26 of sample introduction 
capillary 22 was positioned -20 cm above outlet 23. A vacuum regulator (not 
shown) was used to apply a partial vacuum on the inlet end of the capillary to 
control the sample introduction rate. This configuration resulted in a sample 
stream diameter of approximately 10-pm. Sheath fluid 1§ was ultrapure 
water, delivered to the top of flow cell IS by pump 32, which is a syringe pump 
in the experimental embodiment, but could be any controllable fluid delivery 
system, such as gravity feed, a positive displacement pump or a roller-type 
pump, A second pump 24, which was also a syringe pump in the 
experimental embodiment, was used to move sheath fluid 16 through flow cell 
IS at a volumetric flow rate of nominally 25 - 50 p X min\ 

Fluorescent impurities present in sheath fluid IS were photolyzed by 
passing the sheath fluid through a -1 m-long photolysis cell 2fi before it 
entered the flow cell, as described in U.S. Patent Application 08/727,841 , 
incorporated herein by reference. Fluorescence was induced by irradiating 
the sample stream -25 pm downstream from the capillary tip with a 514.5 nm 
mode-locked argon ion laser beam 2S (82 MHz repetition rate, 200 ps pulse 
width, 20 -30 mW average power), focused to a circular spot of 16 pm (e* 2 ) 
diameter by a 75 mm focal length lens 42. A 60x, 1 2 NA water immersion 
microscope objective 44 was used to collect emitted fluorescence photons 
along an axis orthogonal to the flow and excitation axes. A 600 pm-wide, 
1000 pm -high slit (not shown), with the long axis oriented parallel to the flow 
axis, was placed in the image plane of the collection objective to limit the size 
of the detection volume to ~2 pk. Emitted fluorescence photons were 
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converted to photoelectrons by conventional detection electronics 4fi and 
recorded as photon detection events. 

Hydrodynamic focusing of the sample stream by the sheath fluid to a 
diameter of -1 0 |jm ensured that the entire sample stream flowed through the 
5 central portion of the detection volume. When this is the case, SMD 

efficiencies >90%, essentially limited by the photostability of the fluorescent 
species, are possible. Light passing through the slit was spectrally filtered 
with a bandpass filter (575±15 nm) (not shown) and focused onto a single 
photon counting avalanche photodiode within detection electronics 4S using a 
10 32X microscope objective. Time-correlated single photon counting (TCSPC) 
was used to measure the elapsed time between the excitation laser pulse and 
each detected photon. Background due to Raman and Rayleigh scattering of 
the excitation laser beam by the solvent was suppressed by rejecting photon 
counts detected within -1 ns of the laser pulse. 
15 DATA ANALYSIS 

Each detected photon arriving with a delay >1 ns with respect to the 
excitation laser pulse (gated photon) was recorded. The record consists of 
detection times (with respect to the previous gated photon) for each gated 
photon as well as the arrival time of the photon with respect to the excitation 
20 laser pulse. A photon burst is evidenced by a series of successive gated 
photons recorded at a high rate (40-100 kHz) compared to the background 
counting rate (-5 kHz). Details of the TCSPC apparatus and the algorithm 
used to search data sets for photon bursts is described in, e.g., P.M. Goodwin 
et al., 1895 P. Soc. Photo-opt Ins., pp. 79-89 (1993), incorporated herein by 
25 reference. 

Here, burst search threshold times of 75 and 50 ps were used for data 
collected with 20 and 30 mW of average excitation laser power, respectively. 
Successive gated photons (photoelectrons) detected at time intervals less 
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than the threshold time comprise a photon burst. Each detected photon burst 
is characterized by three parameters: (1) the number of photoelectrons (PE) 
comprising the burst (burst size); (2) the duration of the burst (accumulated 
time below threshold); and (3) the intra-burst fluorescence decay rate (T). The 
5 intra-burst fluorescence decay rate is the reciprocal of the fluorescence 

lifetime r f , where T-\jt f . While the decay rate is the preferred parameter 

for the present invention, it is intended herein that the use of fluorescence 
decay rates include the use of fluorescence lifetimes, if desired by an 
experimenter. The intra-burst fluorescence decay rate was estimated from the 

10 distribution of intra-burst photon arrival times measured with respect to the 

excitation laser pulse, using a maximum likelihood estimator for a background 
free, single exponential decay (J. Tellinghuisen et al. v "Analysis of 
Fluorescence Lifetime Data for Single Rhodamine Molecules in Flowing 
Sample Streams," 66 Anal. Chem., pp. 64-72 (1994), incorporated herein by 

15 reference). To reduce contributions due to background bursts and accidental 
coincidences (a fluorescence burst recorded with two or more analyte 
molecules in the detection volume simultaneously), bursts were time-filtered, 
that is, bursts with durations significantly shorter or longer (see below) than 
the mean molecular transit time across the detection volume were discarded. 

20 The following procedure was used to estimate the probabilities for 

misidentification of the two fluorophores and to characterize the most probable 
sources of error. A Monte Carlo (MC) simulation was used to generate 
fluorescence burst data from a simulated flow cytometry-based SMD 
experiment, given the photophysical properties of the analyte molecules 

25 (absorption cross section; fluorescence quantum yield; fluorescence lifetime; 
photodestruction quantum yield; optical saturation intensity); analyte 
fluorescence emission detection efficiencies of the apparatus; excitation laser 
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parameters (pulse repetition rate, average power, focused spot size); and flow 
parameters (analyte concentration, sample stream flow velocity). The 
simulation accounts for spatial variations of the excitation laser intensity and 
fluorescence collection efficiency within the detection volume, as well as 
diffusion, photobleaching, and optical saturation of the analyte molecules 
during their transit through the detection volume. Good agreement between 
the simulation and the experimental data was achieved by making slight 
adjustments to the estimated optical collection/detection efficiency of the 
apparatus. For the R6G simulation, it was also necessary to adjust the 
photodestruction efficiency, as further discussed in Van Orden et al„ supra. 

P.M. Goodwin, et al., 1895 P. Soc. Photo-op Ins., supra., describes the 
application of a Monte Carlo simulation to model the fluorescence burst size 
distribution detected from single Rhodamine-1 10 molecules excited in flow. 
The simulation has been modified to model the TCSPC data of the simulated 
bursts as well as the fluorescence burst sizes. For each detected photon, a 
random deviate drawn from the appropriate TCSPC arrival time distribution 
(R6G, TRITC, or background fluorescence) was used to assign a time 
separation between that detected photon and the excitation laser pulse. The 
background photon arrival time distribution was compiled directly from the 
experimental data by excluding photons detected inside of photon bursts. 
R6G and TRITC fluorescence photon arrival time distributions were 
constructed by multiplication of single exponential decays (4.0 and 2.3 ns for 
R6G and TRITC, respectively, taken from bulk TCSPC measurements of 
these dyes) by the TCSPC photon arrival distribution measured from an 
uncorrelated light source (ideally, a flat distribution) to account for 
nonlinearities of the TCSPC system. 
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RESULTS AND DISCUSSION 
For all data sets presented in Figures 2A-D, 3A-B, and 4A-F, the 
average excitation laser power was 20 mW. The sheath fluid volumetric flow 
rate was 40 \ik min* 1 ; the e' 2 transit time across the detection volume, derived 
from the autocorrelation function of the photon burst data t was 1 .2 ms; and 
the velocity of the sample stream through the detection volume, estimated 
from the MC simulation, was 1.2 cm/sec. 

Figures 2A-D display 500 ms of raw fluorescence burst data binned at 
250 us intervals from a photolyzed ultrapure water blank (Figure 2A) and 
dilute sample streams of TRITC (Figure 2B), R6G (Figure 2C), and a mixture 
of R6G and TRITC (Figure 2D). A comparison of Figures 2B and 2C clearly 
shows that the fluorescence burst sizes detected for single TRITC molecules 
are, on average, smaller than those detected for R6G. As discussed above, 
the fluorescence burst size of a single fluorophore is determined by such 
photophysical properties as the absorption cross section at the excitation 
wavelength, the fluorescence quantum yield, the fluorescence lifetime, the 
overlap of the fluorescence emissions, and the spectral bandwidth of the 
detection channel. For TRITC, the peak in the absorption spectrum occurs 
near 554 nm. Therefore, excitation at 514.5 nm does not occur as efficiently 
as for R6G, with an absorption maximum at 528 nm. Furthermore, the 
fluorescence quantum yield of R6G (0.9) is almost three times that of TRITC 
(0.35). To improve the relative fluorescence emission detection efficiency for 
TRITC, a spectral bandpass filter was used in the detection channel that 
favored the fluorescence emission maximum of TRITC. 

Figures 3A and 3B display fluorescence burst duration distributions 
(BDD) compiled from data collected from dilute sample streams of TRITC 
(Figure 3A) and R6G (Figure 3B). Experimental conditions were the same as 
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for the data shown in Figures 2A-D. BDDs compiled from data collected with 
the sample stream ON are shown with solid circles; BDDs collected with the 
sample stream OFF and misaligned with respect to the detection volume are 
shown with open circles. The peak in each distribution gives the mean 
detected burst duration for unphotobleached TRITC and R6G molecules 
under these experimental conditions. Note that the mean detected burst 
durations are significantly shorter than the e* 2 transit time (1 .2 ms) derived 
from the autocorrelation function. This is a consequence of the burst search 
algorithm threshold, The majority of single molecule fluorescence bursts fall in 
a time range (0.42-1.38 ms) denoted by the vertical lines and arrows in 
Figures 3A and 3B. Shorter bursts are primarily due to photoelectrons 
associated with Raman scattering that "leak" through the time gates and 
longer bursts are mostly accidental coincidences. Thus, a time-filter can be 
selected to discriminate against background emissions and accidental 
coincidence bursts in the data in real-time, during data collection or, as is 
done here, in a post-processing step. 

Figures 4A-F display size distributions of emitted fluorescence photon 
bursts compiled from data collected under the same conditions as for Figures 
2A-D and 3A-B. Burst size distributions (BSD) shown with closed circles in 
Figures 4A, 4C, and 4E were compiled from data collected from dilute sample 
streams of R6G f TRITC, and a R6G/TRITC mixture, respectively. Open 
circles are BSDs compiled from the water blank (data collected with the 
sample stream OFF and misaligned with respect to the detection volume). 

Figures 4B, 4D, and 4F f display BSDs that were obtained after the 
burst data was time-filtered to remove bursts less than 0.42 ms and greater 
than 1.38 ms in duration. This time-filtering procedure was done to reduce the 
number of background bursts (<0.42 ms) and bursts resulting from two or 
more molecules simultaneously present (>1.38 ms) in the detection volume 



SUBSTITUTE SHEET (RULE 26) 



WO 99/24797 



PCT/US98/24118 



13 

from the data set. Note the marked reduction in the number of smaller (<20 
PE) bursts primarily due to background emissions. It must be emphasized 
that the background BSDs were not subtracted from the analyte BSDs to 
obtain the time-filtered distributions shown in Figures 4B t 4D, and 4F. 

The peaks in the BSDs at -30 and -60 PE correspond to the average 
fluorescence burst sizes detected from single TRITC and R6G molecules, 
respectively, that did not photobleach while crossing the detection volume. 
The fact that peaks are observed in the BSDs indicates that the diameter of 
the sample stream was indeed smaller than the detection volume and was 
well aligned to (he laser beam. Were this not the case, the BSD would simply 
decrease monotonically away from zero PE, i.e., the BSD is a maximum at 0 
PE. 

Figures 4E and 4F show raw and time-filtered BSDs obtained from the 
mixture of TRITC and R6G. Two well resolved peaks at 30 and 60 PE are 
observed, corresponding to different BSDs for TRITC and R6G, respectively. 
This result confirms that these fluorophores can be distinguished at the single 
molecule level based on their different fluorescence burst sizes. The solid 
curves in Figures 4A-4F are BSDs obtained from synthetic data generated by 
the MC simulation described above. The simulation parameters were 
adjusted to give the best fits to the mixture in Figures 4E and 4F. These same 
parameters were then used for the simulations of the pure R6G and TRITC 
solutions. 

Table 1 presents the parameters used in the simulation. Asterisks 
indicate parameters that were adjusted to achieve the best fit. A time-filtered 
burst size detection threshold of 14 PE was chosen as a good compromise 
between background burst rate and SMD efficiency. Time-filtered bursts 
smaller than 14 PE in size are not counted as single molecule fluorescence 
bursts. For this threshold the simulation indicates that 96% of the TRITC 
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molecules and 78% of the R6G molecules that are introduced into the sheath 
flow are detected by apparatus used herein. The smaller efficiency obtained 
for R6G was due to the larger photodestruction quantum yield for this 
molecule. 



Table 1 . Monte Carlo Simulation Parameters 



Parameter 


R6G 


TRITC 


Absorption cross section CT 514nm (10* 16 cm 2 ) 


2.2 


1.6 


riuorescence quaniurn yieiu 


u.y 


U.oO 


Fit inrp^rpnrp lifptimp 






Photodestruction quantum yield (10" 6 ) 


80' 


6 


Saturation intensity (10 5 W cm" 2 s" 1 ) 


0.55 


1.14 


Diffusion constant (10~ 6 cm 2 s 1 ) 


300 


300 


Bandpass filter fluorescence emission transmission 


0.16 


0.27 


Time-gate fluorescence emission transmission 


0.84 


0.74 


Excitation laser e' 2 beam diameter (10^ cm) 


16 


16 


Average excitation laser power (10' 3 W) 


20* , 30" 


20* , 30" 


Sample flow velocity (cm s* 1 ) 


1 .23* , 
0.9" 


1.23*. 
0.9" 


O/erall optical collection/detection efficiency excluding 
bandpass filter and time gate transmissions 


0.093" . 
0.085*" 


0.093" , 
0.085*" 


Detection electronics dead time (10* 6 ) 


4 


4 



* - values adjusted to fit experimental data 

+ - For data in Figures 2-4; ++ - For data in Figure 5 

In order to fit the R6G burst size distribution (Figure 3A) the 



photodestruction quantum efficiency of R6G was adjusted to a final value of 
8x10* 5 . This value is larger than the value of 2x10* 5 reported in the literature. 
In Figure 4B, a shoulder is observed toward the low photoelectron side of the 
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peak. This shoulder is due to R6G molecules that photobleached during their 
transit through the detection volume. A much smaller shoulder is seen for the 
TRITC peak in Figure 4D, due to the higher photostability of TRITC. The 
simulation herein indicates that 44% of the detected R6G molecules 
photobleach while crossing the detection volume compared to 4% of the 
TRITC molecules, 

A single molecule identification (SMI) threshold of 45 PE was chosen; 
bursts greater than or equal to this threshold are identified as R6G single 
molecule fluorescence bursts; bursts <45 PE are identified as TRITC bursts. 
According to the MC simulation, 75% of detected R6G single molecule bursts 
are >45 PE in size and 25% are detected below this threshold. Of the 
detected TRITC single molecule bursts, 96% are <45 PE in size and 4% give 
bursts >45 PE; 97% of the unphotobleached R6G molecules are detected 
with bursts >45 PE. Clearly, photodestruction of R6G molecules as they 
transit the detection volume increases error rates for SMI based on 
fluorescence burst size alone. But, as indicated by Figures 4E and 4F, 
increasing the photostability of the brighter component will reduce the overlap 
of the burst size distributions and decrease the SMI error rates. 

To improve the accuracy of SMI, simultaneous, correlated 
measurements of the burst size and intra-burst fluorescence decay rate were 
performed for each detected burst in R6G/TRITC mixtures. For these 
measurements, the flow velocity of the sample stream through the detection 
volume was lowered to 0.9 cm s' 1 and the average excitation power was 
raised to 30 mW to increase the average number of PE detected per burst. 

Figures 5A and 5B show correlated burst size (PE) and fluorescence 
decay rate (r) measurements from dilute sample streams, containing 
approximately equal amounts of TRITC and R6G (Figure 5A) and containing 
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mainly R6G (Figure 5B). Analyte molecule introduction rates estimated from 
MC simulations are -20 TRITC s* 1 and -20 R6G s' 1 in Figure 5A; and -65 
R6G s" 1 and -3 TRITC s* 1 in Figure 5B. The correlated data are displayed as 
two-dimensional histograms (scatter plots) with darker shades of grey 
indicating increasing numbers of events. Experimental conditions were: a 
sheath volumetric flow rate, 30 \iX min" 1 ; average excitation laser power 30 
mW; sample stream flow velocity, -0.9 cm s°; sample transit time (e~ 2 ), -1 .8 
ms. The burst data were time-filtered (0.25-2.0 ms) to reduce background 
bursts and accidental coincidences; only bursts >20 PE in size are shown. 
Again, no background subtraction was done to arrive at the distributions 
shown in the figures. 

The PE-r plane was divided into four regions (l t ll,lll,IV) for SMI, as 
shown in Figures 5A and 5B. Projections of the correlated data onto the 
fluorescence decay rate and fluorescence burst size axes are shown to the 
left and above each scatter plot, respectively. A and B solid curves plotted on 
the projections are estimates, based on MC simulations of the experiment, of 
TRITC and R6G contributions to the measured distributions, respectively. C 
curves are sums of the estimated TRITC and R6G contributions. Vertical lines 
shown on the projections delimit regions that are used for distinguishing 
TRITC and R6G single molecule fluorescence bursts based on intra-burst 
fluorescence decay rate alone(TRITCsl+IV, r>0.307 ns' 1 ; R6G=il+lll, r<0.307 
ns' 1 ) and burst size (TRITCel+lll, 20<PE<75; R6G=II+IV, PE^75). 

Based on MC simulations for these experimental conditions and a time- 
filtered burst size detection threshold of >20 PE, the estimated SMD 
efficiencies are 97% for TRITC and 80% for R6G. Again, the lower detection 
efficiency for R6G is due to its high photodestruction quantum yield; 
approximately 66% of the R6G molecules photobleach while crossing the 
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excitation laser, whereas only 7% of the TRITC molecules photobleach. The 
increased photobleaching rates (compared to Figures 4A-F) are due to the 
higher average excitation laser power (1.5x) and longer sample transit time 
(1 .4x) used in this experiment. 

According to the MC simulations, for SMI by burst size alone, 99% of 
the detected TRITC molecules give bursts >20 and <75 PE in size and fall into 
Region Mil; only 1% of the detected TRITC molecules give bursts >75 PE 
and fall into Region ll+IV. 53% of the detected R6G molecules give bursts 
>75 PE and fall into Region ll+IV and 47% of the detected R6G molecules 
give bursts >20 and <75 PE in size and fall into Region l+lll. Again, a large 
fraction of the detected R6G molecules are misidentified. This is a 
consequence of detected bursts from photobleached R6G molecules that fall 
into the size range (l+lll) expected for TRITC bursts. 

For SMI by intra-burst fluorescence decay rate (/) alone, 89% of the 
detected TRITC molecules give bursts with r>0.307 ns 1 and 11% give bursts 
with r< 0.307 ns* 1 . For R6G, 88% of the detected molecules give bursts with 
r<0.307 ns* 1 and 12% give bursts with r>0.307 ns"\ Here the errors are a 
consequence of the overlap of the estimated fluorescence decay rate 
distributions for TRITC and R6G. The widths of these distributions are 
determined by the number of PE used to estimate the decay rate. 
Improvement of the overall photon collection/detection efficiency would result 
in more PEs detected per burst and narrower fluorescence decay rate 
distributions, thereby increasing the accuracy of SMI. Photobleaching does 
broaden the R6G distribution indirectly since, on average, fewer PEs are used 
to estimate fluorescence decay rate of photobleached R6G molecules. 
However, this effect is small compared to the direct effect that photobleaching 
has on the R6G burst size distribution. 
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For correlated, two-parameter SMI, two islands corresponding to 
fluorescence bursts detected from unphotobleached TRITC and R6G 
molecules are clearly visible in Figure 5A. TRITC bursts with a mean burst 
size of 53 PE and a mean intra-burst fluorescence decay rate of 0.38 ns* 1 
cluster in Region I; unphotobleached R6G bursts with a mean burst size of 96 
PE and a mean intra-burst fluorescence decay rate of 0.25 ns' 1 cluster in 
Region II. Region III contains mostly bursts from R6G molecules that 
photobleached while crossing the detection volume. This is clear from Figure 
5B where the sample consisted mainly of R6G. Photobleached R6G molecule 
fluorescence bursts comprise a "tail" on the unphotobleached R6G distribution 
that extends from Region II across Region III. 

Based on the MC simulation, probabilities are assigned for detected 
TRITC and R6G bursts to fall into Regions l-IV. These are tabulated in Table 
2. 



Table 2. Single Molecule Identification Probabilities 





1 


II 


III 


IV 


TRITC 


0.883 


0.001 


0.108 


0.008 


R6G 


0.078 


0.489 


0.390 


0.044 



Probabilities for detection in regions delineated for SMI by burst size 
(l+IIUI+IV) and intra-burst fluorescence decay rate (MV,II+III) alone were 
obtained by summing the appropriate entries in Table 2. These calculations 
ignore accidental coincidences and are strictly valid only for SMD rates small 
compared to the reciprocal of the molecular transit time across the detection 
volume. 

It is interesting to calculate the confidence level for identification as 
R6G or TRITC for events that occur in the specific regions. This level 
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depends on the relative detection rates of the two molecules. For example, 
for SMI by burst size alone, the magnitude of the tail of the TRITC BSD 
extending into Region ll+IV depends on the detection rate of TRITC. Likewise 
the magnitude of the shoulder of the R6G BSD into Region I + III depends on 
the detection rate of R6G. 

The fraction of bursts falling into a given region due to TRITC or R6G 
can be calculated from the relative detection rates for these species 
and the SMI probabilities given in Table 2. Table 3 is a tabulation of the 
fraction of bursts detected in a region due to a specific fluorophore for three 
different sample stream concentration ratios of TRITC and R6G 
([TRITC]/[R6G] = 0.25, 1 , 4). For example, comparison of the results for an 
equimolar sample shows that the highest SMI confidences are obtained for 
the correlated measurement in Regions I (93% of the bursts detected in this 
region are due to TRITC) and II (99.8% for R6G). SMI by burst size alone 
gives only a 72% confidence for TRITC (Region I + III) and a 98% confidence 
for R6G (Region II + IV). SMI by intra-burst fluorescence decay rate alone 
gives a 90% identification confidence for TRITC (Region I + IV) and a 87% 
confidence for R6G (Region II + III). 



SUBSTITUTE SHEET (RULE 26) 



WO 99/24797 PCT/US98/24118 



20 

Table 3. Single Molecule Identification Confidence Levels 



1 Parameter - Burst Size 


TRITC|R6G 




kill 


ll+IV 






0.25 


TRITC 


0.39 


0.006 




■ 




R6G 


0.61 


0.994 






1.0 


TRITC 


0.72 


0.02 








R6G 


0.28 


0.98 






4.0 


TRITC 


0.91 


0.08 








R6G 


0.09 


0.92 






1 Parameter: Fluorescence Decay Rate 


TRITC|R6G 




l+IV 


ll+lll 






0.25 


TRITC 


0.69 


0.04 








R6G 


0.31 


0.96 






1.0 


TRITC 


0.90 


0.13 








R6G 


0.10 


0.87 






4.0 


TRITC 


0.97 


0.38 








R6G 


0.03 


0.62 






2 Parameter: Burst Size + Fluorescence Decay Rate 


TRITC|R6G 




I 


II 


III 


IV 


0.25 


TRITC 


0.77 


0.001 


0.08 


0.05 




R6G 


0.23 


0.999 


0.92 


0.95 


1.0 


TRITC 


0.93 


0.002 


0.25 


0.18 




R6G 


0.07 


0.998 


0.75 


0.82 


4.0 


TRITC 


0.98 


0.01 


0.57 


0.47 




R6G 


0.02 


0.99 


0.43 


0.53 
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CONCLUSIONS 

Thus, in accordance with the present invention, single molecules in 
mixed sample streams may be distinguished with high accuracy, based on 
differences in fluorescence burst intensity or a combination of fluorescence 
burst intensity and intra-burst fluorescence decay rate, both using only a 
single excitation wavelength and a single fluorescence emission detection 
channel. For the analytes studied here, TRITC and R6G, SMI confidence 
levels for TRITC are limited mainly by the photostability of R6G. Increasing 
the photostability of the analytes by, for example, the addition of an anti-fade 
reagent is expected to result in significant improvement of SMI confidence 
levels. 

An immediate application is to single molecule DNA sequencing. 
Sequencing multiple identical DNA strands individually and combining the 
results to form a consensus sequence will reduce error rates below those 
reported here. For example, given a random error for sequencing an 
individual DNA strand of 20% per base, the calculations show that the error 
for a consensus sequence having only 10 DNA strands is reduced to less than 
1%. 

While the above examples are based on the pair of dyes Rhodamine 
6G and TRITC, many other dye selections can be identified readily by a 
person skilled in the art of fluorescence detection. Since only one excitation 
wavelength and a single detection channel are used, the various fluorescent 
species chosen for SMI must have overlapping excitation and overlapping . 
emission spectra. The excitation wavelength and emission bandpass filter 
wavelength must be chosen such that the fluorescent species can both be 
detected at the single molecule level with photon bursts of sufficient size 
(greater than about 25 detected photons per burst) for intra-burst fluorescence 
decay measurements to be obtained. 
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In the absence of optical saturation and photobleaching the intensity of 
photon bursts from a given species scales as the product of the absorption 
cross-section and the fluorescence quantum yield of the molecule. 
Photobleaching will reduce the average burst size and broaden the burst size 
distribution. Both the excitation wavelength and detection channel bandpass 
filter center position and width can be adjusted to optimize differences 
between the average detected photon burst intensities of the different 
fluorescent species. 

The fluorescence decay rate of a given species is determined by its 
photophysical properties and interactions between it and its immediate 
environment (e.g., the solvent). The chosen species must have fluorescence 
decay rates sufficiently different such that they can be discriminated at the 
single molecule (photon burst) level. Given the current state of the art for 
overall photon collection/detection efficiencies (1%-5%), the decay rate ratio 
of the chosen species should be at least 1 .5 for SMI based on lifetime alone. 
This ratio is expected to decrease as detection technology improves. 
Differences in the species' average photon burst intensities allows SMI with 
smaller decay rate ratios (closer to unity). 

In general, for single molecule identification by fluorescence burst size 
or by a combination of fluorescent burst size and intra-burst fluorescence 
decay rate, a condition is selected to maximize the resolution between the 
species in the burst size and burst size/fluorescence decay rate plane. SMD 
efficiencies and SMI error rates for a given pair of dyes can be predicted using 
a MC simulation that incorporates photophysical constants of the dyes and 
selected experimental parameters. Some possible fluorescent dye pairs for 
SMI according to the above criteria are provided in Table 4. 
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Table 4 Dye Combinations for SMI 



Dye 1 


abs/em (nm) 


T (ns)' 1 


Dye 2 


abs/em (nm) 


r (ns) 1 


Lisamine 
Knoaamine b 


570/590 


0.475 


Texas Red 


595/615 


0.250 

• 


Rhodamine B 


570/590 


0.556 


Oregon 

A^ ^ A A 

Green 514 


511/530 


0.238 


TRITC 


555/580 


0.475 


BODIPY Fl 


505/513 


0.175 


R6G 


528/555 


0.263 


TRITC 


555/580 


0.475 


Cresyl Violet 
Percniorate 


610/620 


0.313 


BODIPY TR 


589/617 


0.185 


TRITC 


555/580 


0.475 


Oregon 
Green 500 


503/522 


0.244 


Rhodamine B 


570/590 


0.556 


Oregon 
Green 500 


503/522 


0.244 


TRITC 


555/580 


0.475 


Rhodamine 

101 
Perchlorate 


570/590 


0.231 


Rhodamine 
700 


652/678 


0.475 


Rhodamine 
800 


682/716 


0.91 


Cresyl Violet 
Perchlorate 


61 0/620 


0.313 


Rhodamine 
B 


570/590 


0.556 


Cy5 


655/675 


0.665 


Cresyl 
Violet 
Perchlorate 


610/620 


0.313 



The foregoing description of the invention has been presented for 
purposes of illustration and description and is not intended to be exhaustive or 
5 to limit the invention to the precise form disclosed, and obviously many 
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modifications and variations are possible in light of the above teaching. The 
embodiments were chosen and described in order to best explain the 
principles of the invention and Its practical application to thereby enable others 
skilled in the art to best utilize the invention in various embodiments and with 
various modifications as are suited to the particular use contemplated. It is 
intended that the scope of the invention be defined by the claims appended 
hereto. 
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WHAT IS CLAIMED IS: 

1 . A method for identifying single fluorescent molecules in a 
sample stream comprising the steps of: 

forming a sample stream containing a dilute mixture of selected single 
molecule fluorophores, wherein each one of said fluorophores is serially 
ordered in said sample stream; 

passing said serially ordered fluorophores one at a time through a 
detection volume; 

illuminating said sample stream within said detection volume with a 
pulsed laser beam at a single excitation wavelength effective to excite said 
fluorophores; 

detecting emitted photons from each one of said fluorophores to output 
corresponding photon detection events; 

determining from said photon detection events a photon burst size and 
an intra-burst fluorescence decay rate for each one of said fluorophores; and 

correlating said burst size and said decay rate with known 
characteristics for each one of said fluorophores to identify each one of said 
fluorophores. 

2. A method according to Claim 1 , wherein said sample stream has 
a diameter less than a diameter of said laser beam intercepting said sample 
stream. 
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3. A method according to Claim 1 , wherein the step of forming said 
sample stream includes the step of directing a sheath fluid about said sample 
stream to form a focused sample stream within said laser beam. 

4. A method according to Claim 1 , further including the steps of: 
photobleaching said sheath fluid; and 

introducing said dilute mixture of single molecule fluorophores within 
said sheath fluid after said photobleaching. 

5. A method according to Claim 1 , further including the step of 
forming a photon record comprising a number of said photon detection events 
occurring within a time-gated window after a predetermined delay with respect 
to an exciting pulse in said pulsed laser beam and a detection time for each 
one of said photon detection events that is relative to an immediately 
preceding photon detection event. 

6. A method according to Claim 1 f wherein the step of determining 
said size of each said photon burst further includes the steps of: 

time filtering each said photon burst to output filtered burst data 
comprising a number of said photon detection events. 

7. A method according to Claim 1 , wherein the step of determining 
said intra-burst fluorescence decay rate further includes the step of estimating 
said intra-burst fluorescence decay rate from arrival times of said photons 
measured with respect to said excitation laser pulses. 

8. A method according to Claim 1 , wherein said fluorophores are 
selected to maximize the resolution between said fluorophores in burst size 
iand fluorescence decay rate space. 

9. A method for identifying single fluorescent molecules comprising 

the steps of: 
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forming a sample stream containing a dilute mixture of selected single 
molecule fluorophores, wherein each one of said fluorophores is serially 
ordered in said sample stream; 

passing said serially ordered fluorophores one at a time through a 
detection volume; 

illuminating said sample stream within said detection volume with a 
laser beam at a single excitation wavelength effective to excite said 
fluorophores; 

detecting emitted photons from each said fluorophore to output 
corresponding photon detection events; 

determining from said photon detection events a photon burst size for 
each one of said fluorophores; and 

correlating said burst size with known characteristics for each one of 
said fluorophores to identify each one of said fluorophores. 

10. A method according to Claim 9, wherein said sample stream has 
a diameter less than a diameter of said laser beam intercepting said sample 
stream. 

11. A method according to Claim 9, wherein the step of forming said 
sample stream includes the step of directing a sheath fluid about said sample 
stream to form a focused sample stream within said laser beam. 

12. A method according to Claim 9, further including the steps of: 
photobleaching said sheath fluid; and 

introducing said dilute mixture of said single molecule fluorophores 
within said sheath fluid after said photobleaching. 

13. A method according to Claim 9, wherein the step of determining 
said size of said photon bursts further includes the step of: 

time filtering each said photon burst to output filtered burst data 
comprising a number of said photon detection events. 
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14. A method according to Claim 9, wherein said fluorophores are 
selected to maximize the resolution between said fluorophores in burst size 
space. 

15. A method according to Claim 9 f wherein said laser beam is a 
pulsed laser beam. 

16. A method according to Claim 15, further including the step of 
forming a photon record comprising a number of said photon detection events 
occurring within a gated window after a predetermined delay with respect to 
an exciting pulse in said pulsed laser beam and a detection time for each one 

5 of said photon detection events that is relative to an immediately preceding 
photon detection event. 

17. A method according to Claim 1 6, wherein the step of 
determining said size of said photon bursts further includes the step of: 

time filtering each said photon burst to output filtered burst data 
comprising a number of said photon detection events. 

18. A method according to Claim 16, wherein said fluorophores are 
selected to maximize the resolution between said fluorophores in burst size 
space. 
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